Second and higher-order quasi-normal modes in binary black hole mergers 
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Black hole (BH) oscillations known as quasi-normal modes (QNMs) are one of the most important 
gravitational wave (GW) sources. We propose that higher perturbative order of QNMs, generated 
by nonlinear gravitational interaction near the BHs, are detectable and worth searching for in 
observations and simulations of binary BH mergers. We calculate the metric perturbations to 
second-order and explicitly regularize the master equation at the horizon and spatial infinity. We 
find that the second-order QNMs have frequencies twice the first-order ones and the GW amplitude 
is up to ~ 10% that of the first-order one. The QNM frequency would also shift blueward up to 
~ 1%. This provides a new test of general relativity as well as a possible distance indicator. 



I. INTRODUCTION 

Direct detections of gravitational waves (GWs) will be- 
come a reality in the near future with current and future 
detectors, such as LIGO, LISA and DECIGO/BBO % 
GWs can not only provide a test of general relativity but 
also open a new window on the universe. 

One of the most important GW sources are the quasi- 
normal modes (QNMs) of black holes (BHs) [2j. QNMs 
are oscillations of the BH metric perturbations, and they 
are damped by emitting GWs. QNM frequencies are 
complex with the real part representing the oscillation 
frequency and the imaginary part representing the damp- 
ing. By observing the QNM frequencies, one can deter- 
mine the mass and angular momentum of spinning BHs. 

The most promising sources that excite QNMs are bi- 
nary BH mergers. For these events QNMs will be de- 
tected with high signal-to-noise ratio (SNR) (e.g., SNR 
~ 10 5 for - 10 s M Q BH mergers at ~ lGpc by LISA) Q 
since a large fraction of energy [E ~ 1% x M for equal 
mass mergers with a total mass M) goes into GWs. The 
merger rate is also estimated to be large enough 0, H| • 
In numerical relativity there has been recently break- 
throughs for calculating the entire phase of binary BH 
mergers @. The result is that the I = 2, m = 2 QNM is 
dominant, carrying away ~ 1% of the initial rest mass of 
the system pj. 

In this paper we show that the higher-order QNMs 
are also prominent in binary BH mergers and they are 
interesting to search for in the observations and simu- 
lations. Here the "higher-order" is with respect to the 
metric perturbations. Although the dimensionless am- 
plitude of the metric perturbations are negligible small 
when we observe them at detectors, they are relatively 
large near the BH, up to 



^/M ~ {E/M) 1/2 ~ 10% 



(1) 



for QNMs with energy E/M ~ 1% [see Eq. {T5J with 
lu ~ M _1 . ip^ denotes a wave function to be discussed 
later]. Hence the generated second-order perturbations 
would have the amplitude ~ 10% of the first-order ones, 
which may be detectable for high SNR events. 



Higher-order QNMs are essentially analogous to the 
anharmonic oscillations discussed by Landau & Lifshitz 
in "Mechanics" Q. In general, an oscillation with small 



amplitude x is described by an equation, 'x 



2 x = 0, 
acos(ujt + tfi). 
Including the 



which gives the first-order solution x — 
Here a and <p are integration constants, 
second-order term with respect to the amplitude 
the equation becomes 

x + uj 2 x = —ax 2 . 



(2) 



With the right hand side as a source term, we obtain a 
successive solution x = acos(ujt + (/))+ x^ where 
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Here an oscillation arises with a frequency 2to. The im- 
portant point is that the second-order oscillation x^ is 
driven by the first-order one and thereby always exists. 

Our result is that the second-order QNMs would also 
have frequencies twice the first-order ones a/ 2 ) = 2cjW 
and the amplitude up to ~ 10% of the first-order ones. 
Since higher-order QNMs always exist, we may test the 
nonlinearity of general relativity. The purpose of this 
paper is to outline the calculations of the second-order 
QNMs for the Schwarzschild BH and clarify the order-of- 
magnitude estimates. More details will be given in the 
forthcoming papers. We use the units c = G = 1, and 
arbitrarily set M = 1 which we can always recover, if 
necessary. 

Although this paper is the first to study second- 
order QNMs, the second-order analysis is pioneered 
by Tomita [9], and the I = 2, m = case in the 
Schwarzschild spacetime is studied by Gleiser et al. [Tel ]. 
It is also extended to cosmology [ll| and the Kerr 
case 113 ■ 



II. FIRST-ORDER 

Let us consider the metric perturbations to second- 
order, 



5 - a + + /i (2) 



(4) 
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where g^ v is the Schwarzschild metric and the super- 
scripts denote a perturbative order. Expanding the Ein- 
stein's vacuum equation we can obtain basic equations 
order by order [l(| [H[ • 

For the first-order, we use the Regge- Wheeler- Zerilli 
formalism 0, Separating angular variables with 

tensor harmonics of indices (£, to), the equations decouple 
to the even (or polar) part with parity (— 1)^ under a 
transformation (9, (f) — > (it — 9, tt + 4>) and the odd (or 
axial) part with parity (— 1) + . Seven equations for the 
even parity part are reduced to a single Zerilli equation, 
and the other three equations for the odd parity part to 
a single Regge- Wheeler equation. The Zerilli equation is 
given by 



dt 2 



Vz(r) = 1- 



d 2 
or 2 



(5) 



2\ 2A 2 (A + l)r 3 + 6A 2 r 2 + 18Xr + 18 
r 3 (Ar + 3) 2 



where r* = r + 21n(r/2 - 1) and A = (£- l)(£ + 2)/2. By 

Fourier transforming, ipgH(t,r) — J e~ lut 4hnu)( r )d JUJ , the 
Zerilli equation gives a one-dimensional scattering prob- 
lem, which is familiar from quantum mechanics. The 
QNMs are obtained by imposing the boundary condi- 
tions with purely ingoing waves ip 



(i) 

Imtu 



(i) 

trauj 



at the 
at in- 



horizon and purely outgoing waves ip j 
finity. Such boundary conditions are satisfied at discrete 
QNM frequencies wi^ that are complex with the imagi- 
nary part representing the damping. There is an infinite 
number of QNMs for each harmonics (£, to) [l6j |. 

All physical quantities for the first-order even parity 
part can be reconstructed from ipim' If we define 



4P = 



A + l 



K. 



(i) 
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III. SECOND-ORDER 



For the second-order, we also separate angular vari- 

(2) 

ables in terms of tensor harmonics. Instead of ip\', we 
introduce a function 



(2) 
Aim 



2M 



Xr + 3M 



dK, 
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r-2M dt u 



(8) 



in the RW gauge for convenience, which is essentially the 



time-derivative of i\) 



(2) 
tm 



The first-order counterpart 



exactly satisfies x\m = d^£m/®t> an< ^ so tne dimensions 
are 4% - 0{M) and - O{M ) (i = 1, 2). The 
equations for the even parity part are reduced to the 
Zerilli equation with a source term, 



dt 2 



a 2 

dr: 



V z {r) 



X?l(t,r) = S tm (t,r) , (9) 



where the source term Si m is quadratic in ijA 

The most dominant second-order mode is the £ — 4, 
m = 4 even parity one. This is because the dominant 
first-order mode is the £ = 2, m = 2 even parity mode 
and hence Sf m is dominated by the product (£ — 2, m — 
2) x (£ = 2, m = 2) in Eq. ©. This gives the to = 4 
mode that is £ = 4 and even parity. Note that what we 
are considering is the particular solution. Homogeneous 
solutions are just proportional to the first-order ones and 
are therefore trivial. 



IV. REGULARIZATION 

Although it is straightforward to calculate the £ = 4, 
to = 4 source term S44 in terms of ip^ , the raw source 
term does not behave well at infinity and is not suitable 
for calculations. We can find S44 ~ 0(r°) at infinity by 
using the expansion 



in the Regge- Wheeler (RW) gauge according to [I?], Ui 
the GW power is given by 
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Note that in the RW gauge the gauge freedom is com- 
pletely fixed and the physical quantities can be expressed 
by the gauge invariant functions in simple differential 
forms. 

In this paper we focus on the most dominant modes 
in binary BH mergers, i.e^ the £ — 2, to = 2 even parity 
mode for the first-order [6(] and its driving second-order 
mode (the £ — 4, to = 4 even parity mode as shown 
below) . Note that we have to specify m since the degen- 
eracy between to breaks at the second-order. 



where Fj = Fj (T_ ) is some function of T_ = t — r* and 
F'j denotes dFj (T_ ) / cff_ . We need a regularized source 
SZ 9 ~ 0(r- 2 ) at least [i.e., the same as the potential 
Vz ~ 0(r~ 2 )]. At the horizon the source behaves well, 
S44 - 0[(r-2)], with 



4l } =F' H + - a F h + =l(r - 2)F ff + 0[(r - 2) 2 ] , (11) 

where Fff = Fh{T + ) is some function of T + = £ + r*. 

We can regularize the source term by using the regu- 
larized function, 
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which satisfies the Zerilli equation ([9]) with a well- 
behaved source term, S^ 9 ~ 0(r~ 2 ) at infinity and 
S^l 9 ~ 0[(r — 2)] at the horizon. Thus we can remove 
an unphysical gauge-dependent divergence. Note that 
such a regularization is not unique, and for example we 
can replace d/dr with d/dt in Eq. ([12]) . The regular- 



ization is equivalent to adding quadratic terms of the 
first-order gauge invariant function to the second-order 
gauge invariant function, so that it preserves the gauge 
invariance [Isj . 

The explicit form of the regularized source term is 
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where "0 = dip^ /dt and ?// 



dip^ /dr. We can now 
solve the regularized Zerilli equation ([9]). With Fourier 
expansions, this provides a two-point boundary value 
problem with purely ingoing boundary condition at the 
horizon and purely outgoing at infinity. The numerical 
calculation will be presented in a forthcoming paper. 



V. DETECTABILITY 

Without solving Eq. ^ numerically, we can find the 
essential properties of the solutions, i.e., the QNM fre- 
and the order-of-magnitude QNM ampli- 



(2) 

quency w 44 ; 



tude. Since the source term S, 
first-order function ibQ) oc e 



[2)reg 
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is quadratic in the 
*, we have y^ reg 



e 2lw 22„* from Eq. ©. Therefore the second-order QNM 
frequencies are twice the first-order ones, 



W 44n — ZuJ 22n > 



(14) 

which differ from any first-order QNM frequencies w^nn- 
By matching the dimensions in both sides of Eq. ©, we 
also estimate the ordcr-of-magnitude amplitude as 



(2)req 
X44 



,(1) 



(15) 



Then we can derive the second-order QNM energy E^ 2 > 
from the first-order one 2?W. For the QNM waveform, 



= 4 1 2 ) (0)e" 



Eq. ([7]) is given by 
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the first-order energy in 
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where jM ~ 1% for equal -mass mergers [6j. By not- 



9^44 /dt ~ 2w22nV'44 ; j we have a similar ex- 
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pression for the second-order energy as 
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where the second equality uses Eqs. (|T5|) and iflBl. In a 
next paper we will calculate the coefficient as E^ 2 ' jM ~ 

Once we know the QNM frequency to and energy E, 
we can obtain the GW energy spectrum [H, [2(| 

cLE _ 167r 2 £;/ 2 |3a;| 3 

W ~ \uj\ 2 [(2vr/ - 3?w) 2 + (9w) 2 ] 2 ' 



(18) 



and then the characteristic amplitude h c h ar (f) with 
Eq. (5.1) of Flanagan and Hughes [H, 
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where D(z) is the luminosity distance 
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FIG. 1: The characteristic GW amplitudes h c h ar (f) of the 
Schwarzschild black hole (BH) quasi-normal modes (QNMs) 
for two equal-mass binary BH mergers of total mass 10 6 Mq 
at redshift z = 5. The first-order QNM has £ = 2, m = 2 and 
energy E {1) /M = 1%. The second and third-order QNMs 
also appear at twice and three times the first-order frequency, 
respectively, with detectable amplitudes. The dashed line is 
the rms noise amplitude h„(f) = <J 5fSh(f) for the space- 
based detector LISA and ultimate DECIGO. The signal-to- 
noise ratio squared for a randomly oriented source is given by 
(SNR) 2 = Jd(lnf)[h char (f)/h n (f)} 2 . We use a cosmology 
(Q m ,tt A ,h) = (0.27,0.73,0.71). We estimate the third-order 
QNM energy by extrapolating the second-order equation (|17p 
as E (3) /M ~ 15(MQu)) 2 (E^/M) 3 . 



VI. DISCUSSIONS 

Our main results are summarized in Fig. [TJ show- 
ing the characteristic GW amplitudes h c har(f) of the 
Schwarzschild BH QNMs for two equal-mass binary BH 
mergers of total mass 10 6 M Q at redshift z = 5, together 
with the rms noise amplitude h n (f) = y^fShif) f° r 
the space-based detector LISA 0, [2(| and ultimate DE- 
CIGO 4j. The first-order QNM has I = 2, m = 2 and 
energy E^'/M = 1%. We estimate the second-order 
QNM with Eqs. (17 -1115 )1 and the third-order QNM en- 
ergy by extrapolating the second-order equation (fT7| 
as E^/M - f5(A/3w) 2 (^( 1 VAf) 3 . We can see that 
the second and third-order QNMs appear at twice and 
three times the first-order frequency, respectively, with 



detectable amplitudes. 

We can in principle identify the higher-order QNMs 
since their frequencies differ from any first-order ones. 
However the actual identification depends on the SNR 
of the observations [20j or the accuracy of the simula- 
tions, tn simulations with current accuracy we may have 
already mistaken the second-order QNMs for the first- 
order ones. For example 2Mw^ = 0.7473 - 0.1779z is 



close to Mlo 



(i) 

4m0 



0.809f8-0.094f6i hi 



In order to proof that the second-order QNMs actu- 
ally exist, we have to find the second-order QNMs di- 
rectly in the numerical simulations. Such simulations are 
challenging because the mesh size should be less than 
~ 1% X M to resolve ~ I % metric perturbations. Even 
f+f D spherical models for the fully self-gravitating case 
have not found the second-order QNMs [2l|. Simula- 
tions of acoustic black holes may be an alternative for 
this purpose [23]. We also need a mathematically rigor- 
ous definition of second-order QNMs like the first-order 
ones that use the Laplace transformation rather than the 
Fourier transformation 0, H3] • 

Future problems include the Kerr BH case, the odd 
parity mode case, and a more solid third-order formula- 
tion. Since the master equation for Kerr BHs also has 
a source term quadratic in the first-order function [l2l ]. 
we may expect similar results. When BHs have no spin 
before mergers, the final spin is a ~ 0.7 @ and hence 
the Kerr effects may not be so large (as inferred from the 
fact that the QNM frequencies shift by only a small fac- 
tor). The odd parity mode appears when BHs have spin 
before mergers. At the third-order the QNM frequen- 
cies will also shift up to (tp^ /Ad) 2 ~ 1% as suggested 
by the anharmonic oscillations in |8( , probably blueward 
because the GW carries away the BH mass. 

The ratio between the first and higher-order QNM 
amplitude include new information about the total GW 
energy E. Since the observed GW amplitude is h ~ 
(E/M) 1 ' 2 (M/r), this could provide a distance indicator. 
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